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Abstract 
Basing on Mellin-Barnes representations and Miller’s transformation, we present the Gel’fand- 
Kapranov-Zelevinsky (GKZ) hypergeometric systems of the four-loop vacuum Feynman integrals 
with arbitrary masses. Through the GKZ hypergeometric systems, the analytical hypergeometric 
series solutions of the four-loop vacuum Feynman integrals with arbitrary masses can be obtained 
in neighborhoods of origin including infinity. The analytical expressions of the four-loop vacuum 
Feynman integrals can be formulated as a linear combination of the corresponding fundamental 


solution systems in certain convergent region. 
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I. INTRODUCTION 


With the improvement of experimental measurement accuracy at the planned future 
colliders [1-6], Feynman integrals need to be calculated beyond two-loop order. Vacuum 
integrals are the important subsets of Feynman integrals, which constitute a main building 
block in asymptotic expansions of Feynman integrals |7, 8] and are also useful for massless 
theories when a propagator mass is introduced as an intermediate infrared regulator [9]. The 
calculation of multi-loop vacuum integrals is a good breakthrough window in the calculation 
of multi-loop Feynman integrals. In this article, we investigate the analytical calculation for 
the four-loop vacuum integrals with arbitrary masses. 

It’s well known that the completely general one-loop integrals are analytically in the 
time-space dimension D = 4 — 2e [10-13]. At the two-loop level, the vacuum integrals have 
been calculated to polylogarithms or equivalent functions [14-18]. The three-loop vacuum 
integrals are also calculated analytically and numerically in some literatures [19-40]. But, 
very few four-loop vacuum integrals are calculated analytically. The vacuum integrals at the 
four-loop level are only calculated analytically for single-mass-scale [41, 42], equal masses 
[43], and reduction [44, 45]. Recently, Feynman integrals using auxiliary mass flow numerical 
method [46, 47], also can be reduced to vacuum integrals, which can be numerical solved 
by further reduction. In order to further improve the computational efficiency and give 
analytical results completely, it is meaningful to explore new analytical calculating method 
of the four-loop vacuum integrals with arbitrary masses. 

During the past decades, Feynman integrals have been considered as the generalized hy- 
pergeometric functions [48-72]. Considering Feynman integrals as the generalized hypergeo- 
metric functions, one finds that the D—module of a Feynman diagram |64, 73] is isomorphic 
to Gel’fand-Kapranov-Zelevinsky (GKZ) D—module [74-78]. GKZ-hypergeometric system- 
s of Feynman integrals with codimension= 0, 1 are presented in Refs. [79, 80] through 
Lee-Pomeransky parametric representations [81]. To construct canonical series solutions 
with suitable independent variables, one should compute the restricted D-module of GKZ- 
hypergeometric system originating from Lee-Pomeransky representations on corresponding 
hyperplane in the parameter space [82-84]. In our previous work, GKZ hypergeometric sys- 
tems of one- and two-loop Feynman diagrams are also obtained [85-87] from Mellin-Barnes 


representations [68, 69], through Miller’s transformation [88, 89]. There are some recent 
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work in GKZ framework of Feynman integrals [90-108]. 

In our previous work, we have given GKZ hypergeometric systems of the Feynman in- 
tegrals of the two-loop vacuum integral [85] and three-loop vacuum integrals [40]. In this 
article, we derive GKZ hypergeometric systems of the four-loop vacuum integrals with ar- 
bitrary masses, basing on Mellin-Barnes representations and Miller’s transformation. The 
generally strategy for analyzing the four-loop vacuum integrals includes three steps here. 
Firstly, we obtain the Mellin-Barnes representation of the vacuum integral. Secondly, we 
find GKZ hypergeometric system of the vacuum integral via Miller’s transformation. Fi- 
nally, analytical hypergeometric series solutions of the vacuum integral are constructed in 
neighborhoods of origin including infinity. The integration constants, i.e. the combination 
coefficients, are determined from the vacuum integral of an ordinary point or some regular 
singularities. 

Our presentation is organized as following. Through the Mellin-Barnes representation and 
Miller’s transformation, we derive the GKZ hypergeometric system of the four-loop vacuum 
integral with five propagates in Sec. II, six propagates in Sec. III, and seven propagates in 
Sec. IV. And then, we construct the analytical hypergeometric series solutions of the GKZ 
system of the four-loop vacuum integrals in Sec. V. At last, the conclusions are summarized 


in Sec. VI, and some formulates are presented in the appendices. 


Il. GKZ HYPERGEOMETRIC SYSTEM OF THE FOUR-LOOP VACUUM INTE- 
GRAL WITH FIVE PROPAGATES 


The general analytic expression for the Feynman integral of the four-loop vacuum diagram 
with five propagates in Fig. 1 is written as 


U, = (A; a | = 1 


2m)? (q? — m?) (q2 — m2) (g? — m?) 


1 
(qa +% +4 +4)? m?l- m2)’ 


x (1) 
[ 
where D = 4—2e is the number of dimensions in dimensional regularization and A,,, denotes 


D D D D 
the renormalization energy scale, q = (q, 3 %); 4. = amp oat opp an The 


Feynman integral is hard to calculate analytically, if all virtual masses are nonzero. So, one 


can extract the virtual masses from the integral, to facilitate further calculation. Adopting 


M, q 


FIG. 1: Four-loop vacuum diagram with five propagators, which m, denotes the mass of the i-th 


particle and q, denotes the momentum. 


the notation of Refs. [68, 69], the Feynman integral of the four-loop vacuum diagram with 


five propagates can be written as 


( 2 yr" +ico 4 

— RE 2 s; _ 

U, = oe oe as| Imi) I(-s,)rd +s,)| a; (2) 
where s = (s,, +++, 5,), ds =ds,ds,ds,ds, and 
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The integral J, just keeps one mass m,, which can be easily calculated analytically. The 
calculation of the integral J, can be seen in Appendix A. 
And then, the Mellin-Barnes representation of the Feynman integral of the four-loop 


vacuum diagram in Eq. (2) can be written as 


—mé An A? .8-2p tice 4 om? é 
Us= canine me) Le SLUG) TO] 
«(TING -1-s))P4- $+ bs yr6-24+ Ys). (4) 


It is well known that negative integers and zero are simple poles of the function T(z). 
As all s, contours are closed to the right in corresponding complex planes, one finds that 
the analytic expression of the Feynman integral can be written as the linear combination 


of generalized hypergeometric functions. Taking the residue of the pole of I'(—s,), (i = 


1,---, 4), we can derive one linear independent term of the integral: 
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with z, = me (Cases 4). 
We adepi the identity 


r(z-n(1-—z+n)= (T(r -— z) = (-)°r/sin rz, (6) 


which originate from the well-known relation T(z) (1 — z) = m/sin7rz. Then, Eq. (5) can 


be written as 
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with 


T, (a; b| x)= 5 3 3 yA ba aig aig ag a (8) 
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where x = (£,,2,,2,,%,), a = (a,,a,) and b = (b,,b,, b},b,) with 
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Here, we just derive one linear independent term of the integral. We still need to derive 
other linear independent terms of the integral. 


Through the adjacent relations of the coefficient unn 
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the difference-differential operators are written as 
4 
(X0, +a,)T,(a, b |x) = a,T,(a+6,;, b |x), G =1,2), 
i=1 
(9,, +b, - DT (a, b |x) = @, - DZ, b-e,, |x), 
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k 

Here e, € R? (j = 1,2) with e, = (1, 0) and e, = (0, 1), e, = (1, 1), e, € 
R‘ (k = 1,---,4) denotes the row vector whose entry is zero except that the k—th entry is 
LCL 1, 1, 1), 0, = z,n denotes the Euler operators, and Q = 0/0x,,, respectively. 


We can define the auxiliary function 
®.(a, b | x, u, v) = uev? Ti (a, b | x), (13) 


with the intermediate variables u = (u,, u,) = (1, 1), v = (v, v,, Vv, ¥,) = (1, 1, 1, 1). 
Through Miller’s transformation [88, 89], the relations are obtained 

v, ®,(a, b | x, u, v) =a,®,(a, b | x ü; Vig Sal). 

v, 2, (a, b |x, u, v) = (b, — 1)%,(a, b 
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x, u, v), (k= ley); (14) 


which naturally induces the notion of GKZ hypergeometric system. 


In addition, the contiguous relations of Eq. (12) are rewritten as 
4 
Osos + ¥,.)®, (a, b | x, u, v) =a,®,(at+e,,, b | x, u, v), 
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u,u,v,0, ®, (a, b | x, u, v) = 6, (a+ e,, b+e,, |X, u, v). (15) 


k 


The operators above together with Ü,» v, define the Lie algebra of the hypergeometric 
systems |88, 89]. Through the transformation 
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Through Eq. (14), one can finally have the GKZ hypergeometric system for the four-loop 


vacuum integral with five propagates: 
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The row vectors of the matrix A, induce the integer sublattice B which can be used to 
construct the formal solutions in hypergeometric series. 
Defining the combined variables 
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we write the solutions satisfying Eq. (18) as 
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FIG. 2: Four-loop vacuum diagram with six propagators for type A, which m, denotes the mass 


of the i-th particle and q, denotes the momentum. 


Here AT = (a,, +++, a,,) denotes a sequence of complex number such that 
A, -&=B,, (23) 
namely, 
a,+ta,+a,+a,+a,=—-a,, A@+a,+a,+a,+a,, =—a,, 


a,-a@,=b,-1, a,-a,=b,-1, a,-a,=b,-1, a,-—a, =b,—1.(24) 


In the following Sec. V, we will show the analytical hypergeometric series solutions solving 


from the GKZ hypergeometric system in Eq. (18). 


Ill. GKZ HYPERGEOMETRIC SYSTEM OF THE FOUR-LOOP VACUUM IN- 
TEGRALS WITH SIX PROPAGATES 


A. Four-loop vacuum diagram with six propagates for type A 


The four-loop vacuum diagrams with six propagates have two topologies, which one can 
be seen in Fig. 2. The general analytic expression for the Feynman integral of the four-loop 


vacuum diagram with six propagates for type A in Fig. 2 is written as 


E (a2 ja d’q 1 
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Through the Mellin-Barnes transformation, the Feynman integral can be written as 
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Similarly, integrating out q, the Mellin-Barnes representation of the Feynman integral of 


(27) 


the four-loop vacuum diagram with six propagates for type A can be written as 
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Taking the residue of the pole of '(—s,), (i = 1,---,5), one can derive one linear inde- 
pendent term of the vacuum integral: 
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Eq. (6), Eq. (29) can be written as 
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where a = (a,, ---,@,) and b= (0,, ---,0,) with 
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and the coefficient A, is 
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In order to proceed with our analysis, we also define the auxiliary function 
®,,(a, b | x, u, v) = uev? T, (a, b |x), (34) 


with the intermediate variables u = (u,,---,u,;) = (1, 1, 1, 1, 1), v = (u,,---,u%), V = 


= (1, 1, 1, 1, 1, 1). Then one can obtain 
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one derives the GKZ hypergeometric system for the four-loop vacuum diagram with six 


propagates for type A: 
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Here, I is an 11 x 11 unit matrix. 
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Correspondingly the dual matrix A,, of A,, is 


Asa = (Ae i i ee ) ` (39) 


Here, I... is a 5 x 5 unit matrix. The row vectors of the matrix A, a induce the integer 
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sublattice B which can be used to construct the formal solutions in hypergeometric series. 


Defining the combined variables 
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we write the solutions satisfying Eq. (37) as 
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where y = (y,, 7, Ys) Æ = (ai, ++, @,,) denotes a sequence of complex number such 


Aj, a= Ba. (42) 


B. Four-loop vacuum diagram with six propagates for type B 


The general analytic expression for the Feynman integral of the four-loop vacuum diagram 
with six propagates for type B in Fig. 3 can be written as 
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FIG. 3: Four-loop vacuum diagram with six propagators for type B, which m, denotes the mass of 


the i-th particle and q, denotes the momentum. 


Integrating out q, the Mellin-Barnes representation of the Feynman integral of the four-loop 


vacuum diagram with six propagates for type B can be written as 
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where s = (s,, ++, 5,). 
Taking the residue of the pole of '(—s,), (i = 1,---,5), one can derive one linear inde- 
pendent term of the integral: 
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and the coefficient A, is 
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In order to proceed with our analysis, we define the auxiliary function 
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v=e, = (1, 1, 1, 1, 1, 1, 1). Then one can obtain 
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one derives the GKZ hypergeometric system for the four-loop vacuum diagram with six 


propagates for type B: 
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Here, I is a 13 x 13 unit matrix. 
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Correspondingly the dual matrix A,, of A,, is 
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where I... is a 5 x 5 unit matrix. The row vectors of the matrix Aon induce the integer 
sublattice B which can be used to construct the formal solutions in hypergeometric series. 


Defining the combined variables 


47% 12%13% 14 — %8%12%13%15 
=e. ae 
ee Pee, 
— 49% 1271316 = 410% 12%13%17 os Z118 (55) 
Yz = ’ Yy, = ’ Ys = ’ 
ie pe ae oe 2,2, 


we write the solutions satisfying Eq. (52) as 
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where y = (y,, 7, Ys) Æ = (ai, ++, Qia) denotes a sequence of complex number such 
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IV. GKZ HYPERGEOMETRIC SYSTEM OF THE FOUR-LOOP VACUUM IN- 
TEGRALS WITH SEVEN PROPAGATES 


A. Four-loop vacuum diagram with seven propagates for type A 


The general analytic expression for the Feynman integral of the four-loop vacuum diagram 


with seven propagates for type A in Fig. 4 can be written as 
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1 
(GG, +a)? = mle, +)? =m = ns)? =e) 


Integrating out q, the Mellin-Barnes representation of the Feynman integral of the four-loop 
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vacuum diagram with seven propagates for type A can be written as 
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FIG. 4: Four-loop vacuum diagram with seven propagators for type A, which m, denotes the mass 


of the i-th particle and q, denotes the momentum. 
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where s = (s1, ++, Sg) 
Taking the residue of the pole of '(—s,), (i = 1,---,6), one can derive one linear inde- 


pendent term of the integral: 
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and the coefficient A, is 
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[I n| [PO +n)]PG, nE + E nT, + X n) 
In order to proceed with our analysis, we define the auxiliary function 
®_,(a, b | x, u, v) = u®v?-*’T (a, b | x), (64) 
with the intermediate variables u = (u,,---,u,), V = (v, 0), U v e, 
(1, 1, 1, 1, 1, 1, 1). Then one can obtain 
aTa la b | x, u, v) =a,®,,(a, b | x, u, v), (jJ = 1487); 
0,,®,,(a, b |x, u, v) = (b, -1)®,,(a, b |x, u, v), (k=1,+,7). (65) 
Through the transformation 
1 ; 
RE (j =1,:,7), Zīpk = Vko (k=1, Os 
Tı T3 
Zis = » “16 T , 
U U, Uz U,V VgV; U U,Uz U, U0 
T, Ta 
217 = , = , 
U Uy Us Ug U3 VgV; U UU; Ug V, UU. 
e aes EE (66) 
19 UUU UVV, O UUU; 


one derives the GKZ hypergeometric system for the four-loop vacuum diagram with seven 


propagates for type A: 


Aza rar = Brar 
where 
1 1 1 1 0 0 O0 —1 0 0 0 0 -1 -1 
1 1 1 1 0 0 0O 0 —iIi 0 0 0 2b =) 
Ae oe 1 1 0 0 1 10 0 0 -1 0 0 -1 -1 
1 1 0 0 1 10 0 0 0 -1 0 -1 -1 
1 1 0000 1 0 0 0 0 0 -1 -1 
1 100000 0 0 0 0 -1 0 0 
ie Oe ee 
B, T = (—a,, ---, —a,, 6,-1, ---, 6, -—1) 


(67) 


Mdi 


M, qa +49, +q, 


Ms, 43 BE q, 


FIG. 5: Four-loop vacuum diagram with seven propagators for type B, which m, denotes the mass 


of the i-th particle and q, denotes the momentum. 


Here, I is a 14 x 14 unit matrix. 


14x14 


Correspondingly the dual matrix A, a of A,, is 


Aza a ew j lezé ) ’ (69) 


where Iṣẹ is a 6 x 6 unit matrix. The row vectors of the matrix A,, induce the integer 


6 
sublattice B which can be used to construct the formal solutions in hypergeometric series. 


Defining the combined variables 


= 23% 13714715 — 4971371416 — Zioi3 i47 
= ti - , M 
! Z1 ZoZ3Ž4 Zioa ča Z1 Zas Z6 
— 1113214718 — 213% 14% 19 — Z12220 70 
Yı = ’ Ys = , Ye i , ( ) 
Zi Zas Z6 Zior Zi Ža 


we write the solutions satisfying Eq. (67) as 


®,,(z) = (Il 2) 4) (71) 


where y = (y,, «+, Ye), anda? = (a, a,, «++, Qa) denotes a sequence of complex number 


such that 


Ansa Be. (72) 
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B. Four-loop vacuum diagram with seven propagates for type B 


The general analytic expression for the Feynman integral of the four-loop vacuum diagram 


with seven propagates for type B in Fig. 5 is 


Zz D 
7B os: “| : A 2 2 - 2] (Q2 2 
(2r)? (0 — mA + 42)? — m — m) 
1 
(FaFa) = ta a e — m7) 


x (73) 
| 
Integrating out q, the Mellin-Barnes representation of the Feynman integral of the four-loop 


vacuum diagram with seven propagates for type B can be written as 


—m2 2 6 
m An At. 


8-2D pttco m2 s; 2 D 
Ue Ganga (me) Lie SL Ga) Poo] - 1-80) 


. 2 
100 FER m? 


6 


D 3D 4 
3 717 $)M(6- = + ys, )I(7-2D4+ 5° s,)T(4—-D+°s,) 
i=1 i=1 


i=l 


TÈ EE A 


3D a 3 D 
Ne oa Le =a È s (2-5 +s +s, (1+ s, lls, 


i=l i=] 


x 1 z z : (74) 
i=1 i=1 i=1 
where s = (S,, «++, S) 
Taking the residue of the pole of '(—s,), (i = 1,---,6), one can derive one linear inde- 


pendent term of the integral: 


m2 ,4nrA2_\8s-2pD 74 
T: RE 
Os J (4r)8 Ge sin? a2 17a, b | x) ; (75) 
with 
Ta | = Ae, (76) 
n=0 
m2 
where n = (n,, ---, ne), X = (2,, ---, Zoh ©, = za a= (a, ---,a,) and b= (b, ---,8,) 
7 
with 
3D 3D 
a a a, =7— 2D, eae a, =4- D, 
D D 
a,=2- 5, a, =3— D, a, =a, = 1, b, = b, = b, = b, 2-5) 
D 3D 
b = 3-5, b =4—D, b, =5—D, b =6- 5, (77) 
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and the coefficient A, is 


i = (78) 


7B 


®_,(a, b | x, u, v) = u®v? ST (a, b | x), (79) 


with the intermediate variables u = (u,,---,Ug), V = (U,,°°+,Us), U v e 


(1, 1, 1, 1, 1, 1, 1, 1). Then one can obtain 


v, rsa, b | x, u, v) =a,,,(a, b | x, u, v), (j = Lee) 


v, s(a, b |x, u, v) = (b, —1)%,,(a, b |x, u, v), (k=1,---,8). (80) 


Through the transformation 


1 
Zj = ap” (j = 1, ,8), Zs+k Ugo (k = 1, ,8), 
j 
Tı Ta 
Zir 5 > is = ’ 
U Up Uz Uy Us UGV Us UGU Uz U Uy Us UU UGU, U; VgV Vg 
T T 
3 4 
Žig = > %29 = ) 
UUs U3 Uz UU; U6U7 Ug U, Uy U3 U, U3 U7 Ug 
T Xe 
5 6 
Zə = ; » Z2 = ) (81) 
Uy Uy Ug Uz Ug UUV, 


one derives the GKZ hypergeometric system for the four-loop vacuum diagram with seven 


propagates for type B: 


Dog g Vag Dee = Baz ’ (82) 
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1 1 1 1 1 1 0 0 -!1 0 0 0 —1 -l 
1 1 1 1 1 1 0 0 0 —1 0 0 —1 -l 
T 1 1 1 1 0 0O 1 0 0 0 0 0 —1 -l 
A = 
X7B 
1 1 1 1 0 0 0 0 0 0 —1 0 0 0 
1 1 0 0 0 0 0 1 0 0 0 0 0 0 
1 1 0 0 0 0 0 0 0 0 0 —1 0 0 
op 
ve =a Gs , a i 
B,, vce (—a,, > ag; b, 1, ) bs = 1) 
Here, I is a 16 x 16 unit matrix. 


16x16 


Correspondingly the dual matrix Az of A 


where I 


=a 

=i 

= 

—1 i 

=f 

0 
(83) 
(84) 


exs İs a 6 x 6 unit matrix. The row vectors of the matrix A,,, also induce the integer 


sublattice B which can be used to construct the formal solutions in hypergeometric series. 


Defining the combined variables 


y —— Zo Z13 Z14215 1617 pa Z10213214215216718 
1 , 20> 
Zi Bn 23% 485% Z1 9% 3% 4% 5% 
— 413714715 216719 = #1171516 ~20 
Ys = , Yı E 
Zi Zog Z4 Zi Zoa 
= Zis 1621 £ Zi2Ž22 
Ys = ’ Ye => ’ 
Zi Zag 41% 


we write the solutions satisfying Eq. (82) as 


where a? = ( 


Qi, °**, Q,) denotes a sequence of complex number such that 


A, = Bp. 


21 


(87) 


V. THE HYPERGEOMETRIC SERIES SOLUTIONS OF THE FOUR-LOOP VAC- 
UUM INTEGRALS 


A. The general case of the four-loop vacuum integral with five propagates 


In this subsection, we will show the hypergeometric series solutions of the GKZ hyperge- 
ometric system of the four-loop vacuum integral with five propagates with arbitrary masses. 
To construct the hypergeometric series solutions of the GKZ hypergeometric system of the 
four-loop vacuum integral with five propagates in Eq. (18) is equivalent to choose a set of 
the linear independent column vectors of the matrix in Eq. (20) which spans the dual space. 
We denote the submatrix composed of the first, third, fourth and fifth column vectors of the 


dual matrix of Eq. (20) as A,,,,, i.e. 


-1 10 0 
-1 01 0 
Aa = (88) 
-1 001 
-1 0 0 0 
Obviously det A,,,, = 1 Æ 0, and 
B345 = A i A, 
1 1 0 0 0 —1 0 0 0 -1 
00100 -1 1 0 0 -i 
= (89) 
0 0 0 1 0 -1 010 -i 
0 0 00 1 -1 0 0 1 -i 


Taking 4 row vectors of the matrix B,,,, as the basis of integer lattice, one constructs the 


1345 
GKZ hypergeometric series solutions in parameter space through choosing the sets of column 
indices I, C [1,---,10] (i = 1,---,16) which are consistent with the basis of integer lattice 
Bie 


We take the set of column indices J, = [2,6,---, 10], i.e. the implement J, = [1,---, 10] \ 


I, = [1,3,4,5]. The choice on the set of indices implies the exponent numbers a, 


I 
Q 
I 


a, =a, =0. Through Eq. (24), one can have 


4 
a, =a,—4a,, a, = Xb, — a, — 4, a, = 1 — b,, 
i=1 
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a, =1-—b,, a, =1— b, a, => b6,-a,—3. (90) 
Combined with Eq. (9), we can have 
= G7 Oy =O =O => 1, My =I. (91) 


According the basis of integer lattice B the corresponding hypergeometric series solution 


1345? 


with quadruple independent variables is written as 


21 %_ \ "a ( %3%7 \ P2 (Z428 \"3 ( Zso Ņ\ "a 
Pi -ffe a Oe R T ee) 


Ze Zio Ze Zio Ze Zio Ze Zio 
“He Sear OKONO 


where the coefficient is 


[1345 


4 4 
cl) (m) = [Mera +a,+n, )PU+a, — niet +a,+n,) 
i=l 
xT(1 +a, +n +a, +n Ean yn) i (93) 


Using the relation in Eq. (6), one can have 


4 4 
[+04 a 3 n, JI (~ao Eg 2 n,) 
Ciu 0) = = » (94) 
[I aJe + o +n, ECL +a, +n,)P0 +o, H NEC +a, +74) 


where we ignore the constant coefficient term rae ETO, And then, through Eq. (91), the 


corresponding hypergeometric series solution can be written as 


OO 

(De sia Pe Ba Bey (1) 

Oi > y’ Y? Ys Y, X C isas] (n) fisas) ’ 
n=0 


Lym (Yi) "2 (Yo \"3 (YaN ra 

f 5) + ) (95) 
ae a F F a 
with the coefficient is 
D 4 4 
T P DD nV + > n,) 
co) (n) = i=1 i=l , (96) 
C1345] 4 D 
4=1 


Here, the convergent region of the hypergeometric function P in Eq. (95) is 


1< jub al< ludh lul < luh lu] < lul}, (97) 


S a 15 Yz; Yz, Ya) 
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which shows that oT is in neighborhood of regular singularity oo. 
According the basis of integer lattice B,,,., we can also obtain other fifteen hypergeometric 
solutions, which the expressions are collected in Appendix B. The sixteen hypergeometric 


series solutions &” whose convergent region is = can constitute a fundamental solution 


[1345 [1345] 


system. The combination coefficients are determined by the value of the scalar integral of 
an ordinary point or some regular singularities. 

Multiplying one of the row vectors of the matrix B,,,, by -1, the induced integer matrix 
can also be chosen as a basis of the integer lattice space of certain hypergeometric series. 


Taking 4 row vectors of the following matrix as the basis of integer lattice, 


B,.., = diag(—1, 1, 1,1) - Bysis 


134 


SA —1 0 0 0O 1 0 0 Q 1 
0 0 1 0 0 sal 1 0 0 -1 
7 , (98) 
0 0 0 1 O0 eb 0 1 0 -1 
0 0 0 0 1 —1i 0 0 1 -1 


one obtains sixteen hypergeometric series solutions Ë (i = 1,---,16) similarly, which 


[1345] 


the expressions are collected in Appendix C. The convergent region of the hypergeometric 


functions ®® — is 
[1345] 


iiss) = (Was Yor Yar Ya) lM] <1, Il <1, lysl <1, ml <1}, (99) 


which shows that N are in neighborhood of regular singularity 0 and can constitute a 
fundamental solution system. 


Taking 4 row vectors of the following matrix as the basis of integer lattice, 


Bag = diag(1, —1, 1,1) - Bss 


134 


1 1 0 00 -1 0 0 0 -1 
0 0 -1 00 1 —1 0 0 1 
= ; (100) 
00 0 1 0 -1 0 1 0 -1 
00 0 01 -1 0 01 -1 
one also obtains sixteen hypergeometric series solutions Ue (i =1,---, 16): 
o) =o (y,O4,), (101) 


[1345] [1345 
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which the expressions can be obtained by interchanging between y, and y, in os The 


convergent region of the hypergeometric functions oN is 


Bis Z (Yrs Yor Yar W) < lul lul < lal, Wel < luh lul < lul}, (102) 


which shows that P are in neighborhood of regular singularity oo and can constitute a 
fundamental solution system. 


Taking 4 row vectors of the following matrix as the basis of integer lattice, 


1345 


B, = diag(1, 1, —1, 1) - B 


134 


110 0 0-10 0 0 -1 
001 0 0 -1 1 0 0 -1 
= (103) 
0 0 0 -1 0 1 0 -1 0 #1 
000 0 1 -1 0 0 1 -1 
one obtains sixteen hypergeometric series solutions a GS 1,355.16): 
2 E py (Y, > Ya) ’ (104) 


which the expressions can be obtained by interchanging between y, and y, in oe The 


convergent region of the hypergeometric functions oe is 


Fists) = (Uo Yor Yar Y< Yel al< lyh lal < lvl ysl < lul}, (105) 


which shows that os are in neighborhood of regular singularity oo and can constitute a 
fundamental solution system. 


Taking 4 row vectors of the following matrix as the basis of integer lattice, 


Baa E diag(1, 1,1, =1) g Boas 

1 100 0 -1 0 0 0 -!1 

0010 0 -1 10 0 -!1 
= 5 (106) 

0001 0 -1 0 1 0 -1 

0 0 0 0 -1 1 0 0 -1 1 

one obtains sixteen hypergeometric series solutions ee (i =1,---,16): 

0) = 00) gU Oe) » (107) 
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which the expressions can be obtained by interchanging between y, and y, in oi The 


convergent region of the hypergeometric functions os is 


Saag T hs Yor Ys Yal < [ush lul < lyst [Yel < lysh lul <lysl}, (108) 


which shows that oe are in neighborhood of regular singularity oo and can constitute a 


fundamental solution system. 


B. The special case of the four-loop vacuum integral with five propagates 


In order to elucidate how to obtain the analytical expression clearly, we assume the two 
nonzero virtual mass for the four-loop vacuum integral with five propagates. The corre- 
sponding integral of the special case for the four-loop vacuum diagram can be expressed as 
a linear combination of those corresponding functionally independent Gauss functions. 

Through Sec. I, the GKZ hypergeometric system in this special case (m, 4 0, m, Æ 


0, m, =m, = m, = 0) can be simplified as 


As g Vs Pe, = Boos ’ (109) 
where the vector of Euler operators is defined as 

7T 

Us = eat v.) ) (110) 


and the matrix A,, is obtained through deleting the 8th, 9th, and 10th columns of the 


matrix A,: 
1 0 0 0 0 0 1 
0 1 0 0 0 0 +41 
0 0 1 0 0 0 -!1 
A, = (111) 
0 0 0 1 0 0 0 
0 0 0 0 1 0 0 
0 0 0 0 0 1 += 0 
Through Eq. (109), one can have the relations 
a,+a,=-a@,, @+a,=-a,, a,-a,=b,-1, (112) 


with a, = 4— 322, a, = 5— 2D, b, = 2 — 8, and the other a, are zero. The dual matrix A, 
of A, is 
A, =(-1 -1 10 0 0 ie (113) 
The integer sublattice B,, is determined by the dual matrix A,, with B,, = A,,. The in- 
teger sublattice B,, implies that the system of fundamental solutions is composed by two lin- 
ear independent hypergeometric functions. We take the set of column indices J, = {1,--- , 6], 
which implies the exponent numbers a, = 0. And then, the corresponding hypergeometric 
series solution can be written as 
4—32 9p) 
5 (y,)= ,F, M aa EN (114) 
S32 


with y, = z, = Mm? /m?, and „F, is Gauss function: 


a, b X (A)n(b)n 
F, r| = Do (115) 


with (a)n =T (a + n)/T (a). We also can take the set of column indices J, = [1,2,4,---, 7], 
which implies the exponent numbers a, = 0. And then, the another hypergeometric series 


solution can be written as 


D 
ee Dame ee 


5 (y,) = (y,)P?1 F, y,| - (116) 


2 


Here, the convergent region of the hypergeometric functions oe (y,) is |y,| < 1. In the 


region |y,| < 1, the integral correspondingly is a linear combination of two fundamental 


solutions: 


alu) = CPO (y,) + CPO (y). (117) 


[51] ~ [51] [51] ~ [51] 

Multiplying one of the row vectors of the integer matrix B,, by -1, the induced integer 
matrix also can be chosen as a basis of the integer lattice space of certain hypergeometric 
series. And the corresponding system of fundamental solutions for is similarly composed by 
two Gauss functions: 


4-32, 3-D}1 
3D _ , 
o®) (u) = (u) ts oF; : ae 


; (118) 


which the convergent region is |y,| > 1. Correspondingly the integral in the region |y,| > 1 


is a linear combination of two fundamental solutions: 


Salu) = COPO (y) +CD (y). (119) 


[51] [51] [51] [51] 


As M? <m?, m, =m, =m, = 0, 


ee ear | dq 1 
i ii (27)P (q? — m?)q2q?(q, + a + qs + q4)? (q2 — m?) 
SiT eee (120) 
where 
To= (A? a d?q 1 
a a (27)? garage (di + Ge + Is +44)? (03 — m5) 
—m® ,4rA2_\8-2D 3D 
= ——5 RE r4- == Pr G= 2D). 121 
(anys m? ) (4— PG ) (121) 
This result indicates 
—m® ,4nrA2_\8-2D 3D 
1) _ m; M R JESS me 
= (days | a )T(4 z E(- 2D). (122) 
As Mm? > m? and m, =m, =m,=0, I, =1,,, +--+: , where 
hea Oe a d?q 1 
a = (27)P (q? — m?)q2q? (q, +q, + qs + qa)? 
—m® 4TA? ~< 8-2D 3D 
Srei RE r4- r65- 2D). 123 
Gna) e )PG—2D) (123) 
This result indicates 
—m® ,4nrA2_\8-2p 3D 
yo om TA p \8-2 x 
= Tans m ) H= O20): (124) 


Actually, the Mellin-Barnes representation of the Feynman integral in this case can be 


obtained as 


-mê ,AnA2_\s-20 1 fto m?s 
= 5 RE LACET 
a Cod m? ) 2ri In È 3) Ms) 
D 3D 
xI(F -1-s,)P4- = +5,)P(6-2D +5,). (125) 


The residue of simple pole of ['(—s,) provides Cn De (y,), that of simple pole of r(D/2 — 
1—s,) provides COD (y,), that of simple pole of T(4— 22 + s, ) provides Cee (y,), and 
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that of simple pole of [(5 — 2D + s,) provides Coo. (y,), respectively. And the residue of 
the simple pole of r(D/2 — 1 — s,) and T(4 — 22 + s,) can induce 


Q 


=m? ,4rA2_\8-2D 3D D 
39) |) n ms E _ _# 
emeu e) T eebu e (026) 


mie 
A conclusion for the other two nonzero virtual mass for the four-loop vacuum integral 
with five propagates, such as m, # 0, m; 4 0, m, =m, = m, = 0 is analogous to that of 


m, #0, m; £ 0, m, = Mm, =M, = Os 


C. The four-loop vacuum integrals with six or seven propagates 


In our previous works [85, 86], we obtain GKZ hypergeometric systems of some one-loop 
and two-loop Feynman integrals, which show the algorithm and the obvious hypergeometric 
series solutions for the one-loop and two-loop Feynman integrals. Recently, the authors of 
the Refs. [72, 103] also give publicly available computer packages MBConicHulls [72] and 
FeynGKZ [103] to compute Feynman integrals in terms of hypergeometric functions, which 
are meaningful to improve computing efficiency. Through the package FeynGKZ [103], they 
give the examples of some one-loop and two-loop Feynman integrals, which are tested analyt- 
ically by our previous work [85, 86], as well as numerically using the package FIESTA [109]. 

Here, we also evaluate the four-loop vacuum integrals with six or seven propagates using 
FeynGKZ [103], which can be seen in the supplementary material. Note that FeynGKZ 
can’t evaluate the four-loop vacuum integrals with five propagates through Mellin-Barnes 
representations and Miller’s transformation. In the supplementary material, we can see that 
the GKZ hypergeometric systems of four-loop vacuum integrals with six or seven propagates 
are in agree with our above results. Using FeynGKZ, we also give some hypergeometric 
series solutions for the four-loop vacuum integrals with six or seven propagates. The series 
solutions from GKZ hypergeometric systems for special case with the two nonzero virtual 
masses and the three nonzero virtual masses are also showed in the supplementary material, 
and tested numerically using FIESTA [109]. One can see that the computing time using the 
hypergeometric series solutions is about O(10~°) times that using FIESTA, which evaluate 
quickly. 

Note that, except above five topologies, the four-loop vacuum integrals still have five 


topologies, which are one topology with seven propagates, two topologies with eight prop- 
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agates, and two topologies with nine propagates [44]. Due that they have complex math- 
ematical structure, the five topologies of the four-loop vacuum integrals can’t obtain the 
GKZ hypergeometric systems, through Mellin-Barnes representations and Miller’s transfor- 
mation. In next work, we will embed the general four-loop vacuum Feynman integrals into 
the subvarieties of Grassmannian manifold [87], to explore more possibilities of the general 


four-loop vacuum Feynman integrals. 


VI. CONCLUSIONS 


Using Mellin-Barnes representation and Miller’s transformation, we derive GKZ hyper- 
geometric systems of the four-loop vacuum integrals with arbitrary masses. The dimension 
of the GKZ hypergeometric system equals the number of independent dimensionless ratios 
among the virtual mass squared. In the neighborhoods of origin including infinity, we can 
obtain analytical hypergeometric series solutions of the four-loop vacuum integrals through 
GKZ hypergeometric systems. The linear independent hypergeometric series solutions whose 
convergent regions have non-empty intersection can constitute a fundamental solution sys- 
tem in a proper subset of the whole parameter space. In certain convergent region, the 
four-loop vacuum integrals can be formulated as a linear combination of the corresponding 
fundamental solution system. The combination coefficients are determined by the vacuum 
integral at some ordinary points or regular singularities, or the Mellin-Barnes representation 
of the vacuum integral. 

Here, we obtain the analytical hypergeometric solutions of the four-loop vacuum integral- 
s in the neighborhoods of origin including infinity, using GKZ hypergeometric systems on 
general manifold. In order to derive the fundamental solution system in neighborhoods of all 
possible regular singularities, we can embed the vacuum integrals in corresponding Grass- 
mannian manifold [87] through their parametrization. To efficiently derive the fundamental 
solution system, next we will embed the four-loop vacuum integrals into the subvarieties of 


Grassmannian using the a-parametric representation. 
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APPENDIX A: THE CALCULATION OF THE INTEGRAL J, 


For the integral J,, firstly we can integrate out q,: 


I 7” J d” q, 1 
a J (Qa)? (g2): (gq, + Gy +q; +)? 
r(2 +s, +s,) 


= r(1 +s r+ s,) i, ee 


dq 1 
i . Al 
i i: (27)? la? + x(1 — x)(q, + qa + q,)2|? T51 t54 ( ) 


Using the well-known integral 


i PG. DAP ees la (A2) 
OPP FAPT APEE. SA l 
one can have 
oy e i(—)PPT(2 -2 +s, +s,) 1 
ao (Am)PPPA s )E(I + 84) [(q, +g, + q,)?? 72t 
1 
x | darwP/?-2-8 (1 — g) P/s, (A3) 
0 
Through Beta function 
: Š -1 _ P(m)P(n) 
B(m, =f B=) e R A4 
(m,n) = f deia = (AA) 
we can have 
p -IOP TR- zts ts- l- s) l- s) 
n (4n)??T(1+s,)[(1+s,)[(D—-2- s, -—8,) 
$ : . (A5) 


Similarly, we can integrate out q,: 


I 1 
rela kaka ea 
A e i s, +s, Bet (2-1-s,)[(D-—2-s, —8,) 


(47)?/2T(1 + s,)T(2 — P+ s +s DEZ — — $, — 8, —§S,) 
1 
Ta a E (ao) 
And then, we integrate out q,: 
I x dq, 1 
ab ae)? (aes aga)? Porat 
4 
i(-)PPT(4 m 3p + 2 s,)I(2 ls s rÈ? FO SiT 8, — 54) 
= 4 
(4n7)??T(1+s,)T(83 -D+s, +s, +s, rD- 4- x S,) 
1 
— (A7) 
Ae Si 
(43) = 
Last, we integrate out q,: 
ee f d?”q, 1 
(a7 — m2) (a7) = 
4 4 
i 5-3245 s; , 1 | 5-2D+} s 4 
sp S (ra) = aD Dae) Dedy 
Together with Eqs. (A5-A8), one can have 
—1 S; 1 \5-2D+ Si D i 
l = (mye ) (ra) i=1 [IIT -1-9)PG +5) | 
3D 4 4 
= FAAEE (A9) 


APPENDIX B: THE HYPERGEOMETRIC SERIES SOLUTIONS OF THE IN- 
TEGER LATTICE B,,,, 


According the basis of integer lattice B,,,,, one can also construct GKZ hypergeometric 


series solutions below. 
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e J, = [2,5,---,8,10], i.e. the implement J, = [1,---,10] \ Z = [1,3,4,9]. The corre- 


sponding hypergeometric series solution is written as 
2 Daw. D 7 1D = 2 
Da =y? Yf Yr De Cha (n) Jusas) G 
4 p å 
DUE RA 7t tEn) 


cl) m) = = = s (Bl) 
> Ir? D D _D 
NIER +n (2 +n +n )P2- 2 n) 


e J, = [2,4,6,7,9,10], i.e. the implement J, = [1,---,10] \ J, = [1,3,5,8]. The corre- 


sponding hypergeometric series solution is written as 
( ) D D D a ( 
BT me a e Pae Dio 3) 
Dis = YP YS Ye > C1345] (n) fisas) } 
n=0 


r+ Sn)P2-2+ Sn) 
co) (n) = = il (BJ 
[T aE +n EE +n E2- 2 +n n) 


a 


e J, = [2,4,---,7,10], ic. the implement J, = [1,---,10] \ Z = [1,3,8,9]. The corre- 


sponding hypergeometric series solution is written as 


4) _, 2-1, D3 < a 
P 5y? Y = Cas) ORA R 
4 4 
a r(2- 2+ EnB- D+} n) 
N J= 4 l D = p (B3) 
PARE p n Py at tae — 2 T n,)T (2 — DSE n,) 


e I, = [2,3,6,8,9, 10], i.e. the implement J, = [1,---,10] \ Z = [1,4,5,7]. The corre- 


sponding hypergeometric series solution is written as 
B a a D B g < 5 
Dius OP Va W? D Cias 0) Fen > 
4 p á 


a= ž > æ 
nira +n, )T(2- F nT Ea +n) 


e J, = [2,3,5,6,8, 10], i.e. the implement J, = [1,---,10] \ J, = [1,4,7,9]. The corre- 


sponding hypergeometric series solution is written as 


CO 
(6) _, 2-1, D-3 (6) 
D = Y? 4 ek sine ) 
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e I, = [2,3,4,6,9, 10], ie. the implement J, = [1,---,10] \ Z = [1,5,7,8]. The corre- 


sponding hypergeometric series solution is written as 


Da g leo) 
i) =y? tyes cf) (n) fasas , 


[1345 3 4 — [1345] 
4 4 
2-24 Yn )T3B-D+ >¥ n,) 
aa) (n) = i=l i=1 (B6) 
[1345] 4 | D D D D 
[Dn] taOs taO ZENO +n,) 


e J, = |2,3,4,5,6, 10], i.e. the implement J, = [1,---,10] \ J, = [1,7,8,9]. The corre- 


sponding hypergeometric series solution is written as 


OO 
(8) — y P-t (8) 
Pa =Y De is (n) firs y 
n= 


4 4 
N(3-D+ On) r(4—- 224+ ¥n,) 
c® (n) = = = . (B7) 
[1345] 4 D É P Fs 
[fhm )2 Ene- 2 +n- +n- R+) 


e I, = [1,6,---,10], i.e. the implement J, = [1,---,10] \ J, = [2,3,4,5]. The corre- 


sponding hypergeometric series solution is written as 


OO 
(9) _, 2-1, 2-1, 2-1 29 (9) 
Pies =y’ Y? Yz Ye Cisas fisas) » 


PAS a ERD n 
9) (n) = ; i=1 i=1 i (B8) 
pi n! r(2 = 2 +n E(2 +n (R +n (F +n) 


e J,, = [1,5,---,8,10], ie. the implement J,, = [1,---,10] \ Zo = [2,3,4,9]. The 


corresponding hypergeometric series solution is written as 


OO 
(10) aar p38 (10) 
Os g y’ YY Y, 3 Ci sas) WS 245 k 


n=0 
4 
T2-2+ > n)T(3-D+ > n,) 
c00) (n) — i=l i=l l (B9) 
[1345] 4 E B . F 
| Il 7, Jre -5 An(n (S Hn- 5 Hn) 


e J,, = [1,4,6,7,9,10], ie. the implement J,, = [1,---,10] \ Z; = [2,3,5,8]. The 


corresponding hypergeometric series solution is written as 


11 D i BPI p-3 = 11 
ae = y? Yy Y, De Chau ) SF vei) J 
4 
a r2-24+ Yn )T8-D+ ¥ n,) 
11 pan i=1 i=1 
e E . (B10) 


[Iu] +n) +n- F nE +n) 


e J,, = [1,4,---,7,10], ie. the implement J,, = [1,---,10] \ Z} = [2,3,8,9]. The 


corresponding hypergeometric series solution is written as 


(12) 2-1, 3P_4 (12) ( 
Plas a y? Y,? 2i A Frisi ’ 


4 
N3-D+ Yn M4- 224+ Yn) 
09 (a) = L = (B11) 
[Ü nT- 2 Hn (2 +n,)P@-2+n,)P(2-2+n,) 
i=1 


e J,, = [1,3,6,8,9,10], i.e. the implement J,, = [1,---,10] \ Za = [2,4,5,7]. The 


corresponding hypergeometric series solution is written as 


D_ 
ou) = = ye “lye < ae Gea noe ? 
4 
r2-24 Sa)PB-D+¥En) 
con) = = . (B12) 
[1345] | _ D E D D D 
nE S 3 Fa R= E Ea + mag +4) 


e J,, = [1,3,5,6,8,10], i.e. the implement J,, = [1,---,10] \ Z, = [2,4,7,9]. The 


corresponding hypergeometric series solution is written as 


(14) Fly TA G4) 
oo =y? ye aan Tinie , 


e I. = [1,3,4,6,9,10], i.e. the implement J,, = [1,---,10] \ Z; = [2,5,7,8]. The 


corresponding hypergeometric series solution is written as 
a5) _ „2-1, 32-4 a5) ( 
oS Yy Y, z 2 N friisi ’ 
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cQ?) (n) — i=1 al À (B14) 


e J,, = [1,3,4,5,6,10], i.e. the implement J,, = [1,---,10] \ Ze = [2,7,8,9]. The 


corresponding hypergeometric series solution is written as 


(16) _ , 2p-5 $ (16) 
ee =y, = Ce Ta eat ’ 
hj Z 
r(4— 224+ Yn) -2D+ Y n,) 
(16) (n) = i=1 i=l , (B15) 


[1345] 4 
1 _D 
fa n,'T(2 2 + n,) 


APPENDIX C: THE HYPERGEOMETRIC SERIES SOLUTIONS OF THE IN- 
TEGER LATTICE B, 


According the basis of integer lattice B. 


i345) ONE Can construct GKZ hypergeometric series 


solutions below. 
e J, = [1,2,7,---,10], i.e. the implement J, = [1,---,10] \ J, = [8,4,5,6]. The corre- 


sponding hypergeometric series solution is written as 


1 D Ije fe D ay = 1 
pl ) = ye y2 y? YP > Coe (0) esas 5 


[1345] 


A E a ae (Cl) 


where the coefficient is 


co) (n) = = i=l l (C2) 
IP +n,) 
=1 


e J, = [1,2,6,---,9], ic. the implement J, = [1,---,10] \ J, = [8,4,5,10]. The corre- 


sponding hypergeometric series solution is written as 


[1345] E 


DY ne 0 E bis 2 
ol ) y? Ye y? 2 (mn) P isa 3 


4 4 
r+ X n) -24 En,) 


e I, = {1,2,5,7,8, 10], i.e. the implement J, = [1,---,10] \ J, = [3,4,6,9]. The corre- 


sponding hypergeometric series solution is written as 


D1 D1 Bay 
g6) =Y? $ A “y? : 2, A (1) fisas ? 


[1345] 


8) (n) = — i=l i=l | (CA) 
T EO +n E +n, E(B +n- 2 +n,) 


e I, = [1,2,5,---,8], ic. the implement J, = [1,---,10] \ J, = [3, 4,9,10]. The corre- 


sponding hypergeometric series solution is written as 


[1345] 


pl AS YP YF 2 Coe (1) fisas i 


e J, = |1,2,4,7,9, 10], i.e. the implement J, = [1,---,10] \ J, = [8,5,6,8]. The corre- 
sponding hypergeometric series solution is written as 


5) =op Pop B] Di 5 
ol ) y? y2 y? Deu (1) fisas 5 


[1345] a 


4 4 
Pe BLOF ee a) 


6) (n) = =r i=1 i (C6) 
n IER +n, E(B +n- 2 +n )0(8 +n) 


e I, = |1,2,4,6,7,9], i.e. the implement J, = [1,---,10] \ J, = [3, 5,8,10]. The corre- 


sponding hypergeometric series solution is written as 


g6) 


[1345] 


D 
2 


OO 
= es REI (6) 
=y Yy 2 C isat (n) Ts X 
n= 


-Bgn BDF 
2 i & Th 


(8) (n) = i i=1 i=1 f (C7) 
[E nO- F HEE + n2- +n +n) 


e J, = |1,2,4,5,7, 10], i.e. the implement J, = [1,---,10] \ Z, = [8,6,8,9]. The corre- 
sponding hypergeometric series solution is written as 


Dı DB 
2 2 


(7) 2-1, P17 
® =y Y, Do “ss (1) fisas) 4 


[1345] 
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e J, = [1,2,4,5,6, 7], ie. the implement J, = [1,---,10] \ J, = [8,8,9, 10]. The corre- 


sponding hypergeometric series solution is written as 


OO 
D 
ps = y?! (8) 7 
a Y, 2L Cisa P fias i 
n= 


4 
N(3-D+ dn, )r(4—- 22+ ¥n,) 
eas ta = = (C9) 
[In] 2 +n )rG En- E +n JQ -F +n) 


e J, = [1,2,3,8,9,10], ie. the implement J, = [1,---,10] \ J, = [4,5,6,7]. The corre- 


sponding hypergeometric series solution is written as 


OO 

(9) _,,@-1,,2-1,,2-1 (9) 

Toy Y, Ys Y, Do “ise (n) fisas , 
n= 


c0) (n) = = = . (C10) 
n IE(B +n,)P2- 2 +n,)0(2 +n,)P(2 +n) 


e J,, = [1,2,3,6,8,9], ie. the implement J,, = [1,---,10] \ J,, = [4,5,7,10]. The 


corresponding hypergeometric series solution is written as 


(10) __ aly =] 
oa Y? yet eA Frags 


2. 4 2 4 
T2-8+ Yn )P-D+ ¥n,) 


c00 (n) = 


[1345] l (C11) 


4 

I n,!|P2@- 2+n,)P2- 2 Hn EB +n,)0(2 +n,) 

e J,, = [1,2,3,5,8,10], ie. the implement J,, = [1,---,10] \ Za = [4,6,7,9]. The 
corresponding hypergeometric series solution is written as 


a1) — 2—1 
rae y? E va en Fisted 


c01) M i=1 =" 
Cr san) 2) = i 5 (C12) 


--,10] \ Z, = [4,7,9,10]. The 


e J,, = [1,2,3,5,6,8], ie. the implement J,, = [1,- 


corresponding hypergeometric series solution is written as 


E —1 
ot one =y? 2 et Í ase ’ 


4 4 
N(3—-D+ Yn )r(4— 22+ ¥ n,) 
(2) (n) = i=l i=1 (C13) 
[1345] 4 D D z a 
[Mag re- F+nr@— 2 A + n)@- 3 +14) 
the implement J,, = [1,---,10] \ 4, = [5,6,7,8]. The 


e I, = [1,2,3,4,9, 10], ie. 
corresponding hypergeometric series solution is written as 


(13) __ =i 
? > y? y 2 3 ae Pisas J 


[1345] 
4 4 
(2-2 +E nB -D+ Yn) 
OEE 5 7 = 
[ TI n! E(R +n E2- 2 +n,)P2@- R nE +n,) 
i=1 
the implement J,, = [1,---,10] \ 4, = [5,7,8,10]. The 


e I, = [1,2,3,4,6,9], i.e. 
corresponding hypergeometric series solution is written as 


= —1 
a ie =y? yet ot Fis 2 


4 
r(3-D+E n4- 2+ Yn) 
c08 (n) ; i=1 i (C15) 
[1345] : = F = F 
pi nire Son ROS 2a Fn +n) 
the implement J,, = [1,---,10] \ Ze = [6,7,8,9]. The 


e I. = [1,2,3,4,5, 10], ie. 
corresponding hypergeometric series solution is written as 


— —1 
a05 = yz D A Jia ’ 


4 
N(3—-D+ ¥n,)r(4— 22+ ¥ n,) 
aa) ~ 74 D = D = D (C16) 
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the implement Je = [1,---,10] \ Ze = [7,8,9,10]. The 


e I; = [1,2,3,4,5,6], i.e. 
corresponding hypergeometric series solution is written as 


go) — Z c16) (1) fisis 


[1345] [1345] 
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